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ABSTRACT 

The 75'0-S'VF5^| 2.5-45 fim infrared spectroscopic observations of the nucleus of the 
Seyfert 2 galaxy NGC 1068 (see companion paper) are combined with a compilation 
of UV to IR narrow emission line data to determine the spectral energy distribution 
(SED) of the obscured extreme-UV continuum that photoionizes the narrow line 
emitting gas in the active galactic nucleus. We search a large grid of gas cloud models 
and SEDs for the combination that best reproduces the observed line fluxes and NLR 
geometry. Our best fit model reproduces the observed line fluxes to better than a 
factor of 2 on average and is in general agreement with the observed NLR geometry. 
It has two gas components that are consistent with a clumpy distribution of dense 
outflowing gas in the center and a more extended distribution of less dense and more 
clumpy gas farther out that has no net outflow. The best fit SED has a deep trough at 
~4 Ryd, which is consistent with an intrinsic Big Blue Bump that is partially absorbed 
by ~ 6 x 10 19 cm -2 of neutral hydrogen interior to the NLR. 



1. Introduction 

The intrinsic spectral energy distribution (SED) of active galactic nuclei (AGN), which 
extends from the radio up to 7-rays, cannot be directly observed from the Lyman limit and up to 
several hundred eV due to Galactic and intrinsic absorption. However, the extreme-UV (EUV) 
and soft X-ray continuum can be investigated indirectly by the infrared coronal line emission. 
These lines are emitted by collisionally excited forbidden fine-structure transitions of highly 
ionized atoms, whose ionization potentials extend well beyond the Lyman limit up to hundreds of 
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eV. Unlike the strong permitted lines of these ions, which are also emitted in the obscured EUV, 
the reddening-insensitive forbidden IR coronal lines and semi-forbidden optical coronal lines can 
be observed. Therefore, when photoionization is the main ionization mechanism, the coronal lines 
can provide information on the intrinsic obscured SED and the accretion process that powers the 
AGN. This information can be extracted by photoionization models of the NLR. 

The coronal lines are collisionally suppressed in the dense broad line region (BLR) close to 
the continuum source and are efficiently emitted only from the more rarefied gas in the narrow 
line region (NLR), hundreds of pc away from the center. It is well established that large quantities 
of gas attenuate the continuum emission in many AGN. These gas clouds are detected by narrow 
UV absorption lines (e.g. Kriss et al. 1992a| ) or by X-ray absorption features and emission lines 



(e.g. George et al. |1998| ). Although their exact location along the line of sight is unknown, there 
are reasons to believe that in some cases they may be inside the NLR. In particular, the warm 
absorbers that block the X-ray continuum appear to cover a large fraction of the continuum 
source (George et al. 1998j ). This raises the possibility that in some AGN the ionizing SED, 



which is traced by the coronal lines, is not the intrinsic one produced by the accretion process, 
but rather one that is filtered by intervening absorbers inside the NLR. It has been proposed that 
such absorbers are common in Seyfert galaxies, and are responsible for the observed correlations 
between the soft X-ray slope and the narrow emission line spectra of Seyfert 1.5 galaxies (Kraemer, 
Ruiz & Crenshaw |1999| ). 



This study of the Seyfert 2 galaxy NGC 1068 is part of the ISO-SWS program on bright 
galactic nuclei. Previous studies in this program include the reconstruction of the SED of the 
Seyfert 2 Circinus galaxy (Moorwood et al. 1996| ; Alexander et al. 1999| ) and of the Seyfert 1 



Galaxy NGC 4151 (Alexander et al. |1999| ). In both cases we found evidence of a "Big Blue Bump" 



signature of a thin accretion disk (Shakura & Sunyaev 1973). However, in the case of NGC 4151 
this structure is masked by a deep absorption trough of an absorber situated between the BLR 
and the NLR, which filters the light that photoionizes the NLR. 

In this paper we apply our SED reconstruction method to NGC 1068, one of the closest, 
brightest and most extensively studied Seyfert 2 galaxies, which is considered a prototype of this 
AGN class. The first detection of broad permitted emission lines in the polarized light of NGC 1068 



(Antonucci & Miller 1985 ) provided a major argument for the Seyfert 1 and 2 unification scheme 
(Antonucci |1993 ). This scheme postulates that the two Seyfert types have both broad and narrow 
line regions and an obscuring torus that lies between the two. When the torus is face on and 
the BLR is directly observed, the AGN is classified as a Seyfert 1. When the BLR is obscured 
by the torus, the AGN is classified as a Seyfert 2, and the BLR can be observed only indirectly 
in scattered polarized light. The factors that determine the accretion properties and the ionizing 
SED of AGN are currently unknown. However, the Seyfert unification picture implies that all 
possibly relevant factors being equal, such as luminosity, host galaxy type or redshift, the intrinsic 
SED of both AGN types should be similar. It is therefore of interest to complement our previous 
study of the nearby luminous Seyfert 1 galaxy NGC 4151 with a corresponding study of a nearby 
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luminous Seyfert 2 galaxy with a similar host galaxy type, such as NGC 1068. The ISO-SWS 



observations of NGC 1068 are presented in a companion paper (Lutz et al. 2000) and are used 



there to derive the gas density and to place constraints on the structure and dynamics of the NLR. 

This paper is organized as follows. In §|| we summarize the physical properties of the nucleus 
of NGC 1068 that are needed for constructing the photoionization models and constraining their 
results. In §|3| we present the emission line flux compilation that we use in our modeling. In §|] we 
briefly discuss the construction and fitting of the NLR photoionization models. We present the 
results in §|5| and discuss them in SS. 



2. The physical properties of NGC 1068 

NGC 1068 is a barred spiral galaxy at z = 0.0036 (distance D = 16.6 Mpc for 
Hq = 65kms _1 Mpc -1 ) with magnitude = 9.17 (e.g. Lipovetsky, Neizvestny & Neizvestnaya 



1988 ). Observations of the nucleus of NGC 1068 and models of the nuclear line emission 
indicate that the gas in the nucleus forms a complex system, which is composed of various 
spatial and dynamical components. These are excited by several physical mechanisms, including 
photoionization by the nuclear continuum, photoionization by hot stars, and possibly also by 
shocks and energetic particles from a radio jet. In order to isolate the effects of the nuclear 
continuum and to construct photoionization models of the NLR it is necessary to understand the 
morphology and content of the galactic nucleus. We present here a brief overview of the properties 
of the nucleus that are relevant to this work. 



2.1. The galactic nucleus 



The most prominent morphological feature in the nucleus of NGC 1068 is the asymmetric 
bi-polar pattern of both the radio and the optical line emission. The radio emission extends 
over ~ 15" and has a sharply defined northern lobe and a weaker diffuse southern lobe (Wilson 



& Ulvestad 1983 ; Muxlow et al. 1996| ). The southern lobe is both smaller and redder, which 
is consistent with the picture that the large northern lobe is observed above the galactic disk, 
generally facing the observer, and the southern lobe is seen through the galactic disk (Unger et al. 



1992 ; Macchetto et al. 1994 ). Images of the NLR in low excitation lines show mainly the northern 
cone (Cecil, Bland & Tully [1978 ; Unger et al. 1992| ; Evans et al. 1991| ). The precise value of 
the opening angle associated with the emission maps depends on the way the edge of the cone is 
defined and on the assumed position of the nucleus. The location of the nucleus can be determined 
to within ~0.05" by the center of symmetry of the UV polarization pattern (Capetti et al. 1995, 



Kishimoto 1999[ ). The position of the nucleus does not appear to coincide with the maximum of 
the continuum emission, which implies that the nucleus is heavily obscured even in the infrared. 
This is consistent with an obscuring torus of column density in excess of 10 24 cm -2 , as is inferred 
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from X-ray (Marshall et al. 1993] ) and CO observations (Tacconi et al. 1994]). The positioning of 
the nucleus makes it possible to estimate the opening angle of the radiation cone at ^ 70° up to 
~100°. The ionization cone appears to be only partially filled. 



Marconi et al. ( |1996| ) find that the coronal line emission peaks 0.5" NE of the nucleus and 
extends up to ~4". They note that the blueshift of the emission line profile centroid increases, 
and the FWHM of the profile decreases with E\ on (the ionization energy required to produce 
the emitting ion from the preceding ionization stage). They interpret this as evidence that the 
emission lines are emitted from outflowing material. The high ionization lines are emitted in the 
inner light cone, where the velocity field is relatively coherent, while the lower ionization lines are 
emitted from slower, more extended areas with different velocities. The NLR [O in] A5007 emission 
extends over the few inner arcseconds (Evans et al. |1991 ; linger et al. 1992| ; Dietrich & Wagner 
1998[ ) . The extended emission line region (EELR) [O in] A5007 emission extends over more than 
10" (Unger et al. |l99^ ). 



The unresolved BLR is seen only in scattered polarized light and has FWHM of ~ 3000 
km s" 1 (Miller et al |l990l) . As is seen in other Seyfert 2 galaxies (Wilson |1988| ; Capetti et al. 
|1996P , the morphology of the NLR emission maps is correlated with that of the radio structure 



(Wilson & Ulvestad 1983; Capetti et al. 1997). In particular, a 3-dimensional reconstruction of 
the positions of individual clouds based on polarimetric measurements also indicates that the 
northern NLR cone is directed towards the observer, and the southern part away from the observer 
(Kishimoto |1999| ). The close correspondence between the NLR and the jet suggests that the jet 
outflow sweeps and compresses the ambient gas and thereby increases the line emissivity. The 



NLR has a very complex structure (Macchetto et al. 1994) and displays large scale dumpiness 



It is composed of many line emitting cloud complexes (Alloin et al. 1983; Meaburn & Pedlar 
Evans et al. 1991; Dietrich &: Wagner |1998 ). The individual components have FWHM 



1986 



ranging from ~200 km s 1 to ~700 km s , and extend over ~2500 km s 1 in velocity space, 
resulting in an integrated [Oiii] A5007 profile with FWHM of 1150 km s _1 (Dietrich &; Wagner 



1998). The velocity field of the NLR clouds with the lowest FWHM is consistent with rotation 
around the nucleus. The bulk velocities of intermediate FWHM clouds are clearly split relative to 
the symmetry axis of the radio jet, which suggests that they are associated with the interaction 
between the radio jet and the NLR gas. The highest FWHM clouds are associated with highly 
polarized emitting structures (Capetti et al |1995|) , and could therefore be a reflected image of an 
inner obscured region. The EELR appears to follow the galactic rotation (Unger et al. 1992j ). 



The large ISO apertures (14" x 20" to 20" x 33") includes both the outflowing inner NLR and 
the rotating EELR. 
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2.2. The line emitting gas 



There is evidence that at least three different ionization mechanisms are at work in the NLR. 
The very high ionization states are probably due to the the central continuum source. Marconi et 
al. ( |1996| ) find that the infrared coronal line ratios point to photoionization as the main excitation 
mechanism of the coronal gas, and that both collisional excitation or photoionization by very hot 
stars can be ruled out. At lower ionization states, the jet / ISM interaction can provide internal 
sources of excitation in addition to the external central continuum, for example by fast shocks 
or cosmic rays. The morphological connection between the radio and line emission suggests that 
the jet outflow shapes the NLR. Estimates of high gas temperatures (Kriss et al. 1992b| ) and the 
existence of dense but highly ionized clouds near the nucleus on both sides of the jet axis suggest 



that the jet may also play a role in photoionizing the clouds (Capetti et al 1997; Axon et al. 1998 



Dietrich & Wagner 1998). Hot stars are a third ionization mechanism. Unresolved UV continuum 
point sources in the inner 7" x 7", which are not observed in [O in] A5007, could be OB associations 



(Macchetto et al. 1994 ). Hot stars are certainly a component in the ring-like structure that 
surrounds the nucleus. An ellipse of Hn emission delineates the NLR at an average radius of 13" 



(Cecil, Bland & Tully 1978 ; Bruhweiler et al. 1991 ) and starburst knots and CO emission encircle 
the nucleus and the NLR at an average radius of 18" (Planesas, Scoville & Myers 1991| ). However, 
the overall similarity in the profiles of the mid-infrared high-excitation lines, which cannot be 
excited by hot stars, and the profiles of the intermediate excitation lines (-Eion > 30 eV), which 
could be excited by hot stars, strongly suggests that gas excited by hot stars within the large 
ISO aperture does not contribute more than ~20% to these NLR line fluxes (Lutz et al. 2000). 



An additional complication due to stars is contamination of the observed emission line spectrum 
with stellar absorption features. This appears in the difference spectrum of the 3O"0 and 18" 
HUT apertures, which shows both a reddened early- type stellar continuum and stellar absorption 
features (Kriss et al. 1992b|) . 

As is discussed by Lutz et al. (pOOOj ) , the ISO-SWS line ratios indicate that the mid-IR lines 
are emitted from gas with a hydrogen density of ~ 2000 cm -3 . The NLR appears to contain also 
higher density gas. The density of individual knots in the high ionization core is estimated at 10 4 
to 3 x 10 cm~ 3 from HST measurements of the [Sn] doublet, while the overall density ~1" from 



the nucleus varies between 10 3 to 4 x 10 3 cm 3 (Capetti et al. 1997j ). 



The ionization parameter^] also appears to vary across the NLR. The 
[Om] A5007/(Ha + [Nil] A6584) ratio (Capetti et al. |1997|) traces a high excitation core 
(log U ~ —2.5) in the inner 1" x 2" north of the nucleus, followed by a lowered ionization 
halo (log U ~ —3.3) out to ~ 4" from the nucleus, and then intermediate ionization filaments 



5 17 = Qion/47rr 2 nc, where Qi on is the ionizing photon luminosity, n is the hydrogen density at the illuminated face 
of the cloud, r is the distance of the face of the cloud from the continuum source and c is the speed of light. The 
local ionization parameter in the cloud falls with increasing depth due to absorption and geometrical dilution of the 
radiation field. 
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(log 17 



-2.8) which extend out to the EELR. 



Estimates of dust reddening in the NLR of NGC 1068 range from -Eb-v = 0.07 for the 
continuum (Kriss et al |1992b|) , to 0.20 (Marconi et al |i"996l) , 0.40 (Shields & Oke [l97|; Neugebauer 
et al ^980] Ward |i~987l ) , and £ B -v = 0.52 (Koski |i~978l) for the NLR. The analysis of Kraemer et al. 



(1998) indicates that there may be some dust mixed with gas in varying amounts. The observed 



line ratios in the NLR suggest that the O/N abundance is less than solar. Netzer (1997) and Netzer 
& Turner ( 1997| ) interpret this as an indication of under-abundant oxygen (see also Sternberg, 
Genzel & Tacconi 1994j ), and propose that the abundances of He:C:N:0:Ne:Mg:Si:S:Ar:Fe 
relative to hydrogen are (100 : 3.7 : 1.1 : 2.7 : 1.1 : 0.37 : 0.35 : 0.16 : 0.037 : 0.4) x 10~ 4 , 
respectively. Kraemer et al. ( |1998j ) interpret the line ratios as showing an overabundance of 
nitrogen, as well as hinting at higher than solar iron and neon, and propose abundance ratios of 
(100 : 3.4 : 3.6 : 6.8 : 2.2 : 0.33 : 0.31 : 0.15 : 0.037 : 0.8) x 10~ 4 . 

Finally, Kraemer et al. ( |1998| ) model the NLR emission with a multi-component gas model, 
and suggest that the low ionization emission lines are emitted by gas that is partially screened 
by a dense, optically thin component. They use their models to estimate that the filling factor is 
F ~ 10" 4 . 



2.3. The ionizing continuum 

Unlike the situation in Seyfert 1 galaxies, where it is possible to observe the intrinsic SED 
outside the obscured range, the intrinsic SED of NGC 1068 cannot be directly observed even in 
the optical or X-ray bands. Only a small fraction of the AGN light is scattered into the line of 



sight, and can be observed against the host galaxy in polarized light. Pier et al ( 1994 ) compiled 
various continuum measurements in the optical, UV and X-ray, and carefully took account of 
aperture differences, star-light contamination, reflection by dust and bremsstrahlung emission 
from the scattering plasma. They conclude that the resulting reflected SED is broadly similar to 
that observed in Seyfert 1 galaxies. Pier et al ( |1994| ) also list various estimates of the fraction of 
light reflected by the scatterer. These values range from f Te Q = 10~ 3 (Bland-Hawthorn & Voit 
1993D to ~0.05 (Bland-Hawthorn, Sokoloski & Cecil |1991| ). Pier et al ( |1994| ) argue that the most 



reliable estimate is / re fl ~ 0.01 to within a factor of a few. Because the reflectors are more than a 
hundred light years away from the continuum source (Miller et al. 199C| ), short-term continuum 



variability is unlikely to affect the reconstruction of the reflected SED. 

We adopt the Pier et al. nuclear continuum SED in the UV and X-ray as the template SED 
(Fig. Q), and further extend it from 10 keV to 100 keV with a slope of F v oc i/ — 1 . We enumerate 
on the unobserved UV to soft X-ray range to find the best fitting SED. 
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3. The observed line flux compilation 



In addition to the ISO-SWS mid-IR lines fluxes presented in Lutz et al. ( pOOC ), we compiled 



a list of UV to IR lines from the literature. The compilation initially included a list of ~ 120 
measured line fluxes, which were obtained over the last ~30 years using various instruments with 
different apertures, spectral resolutions and reduction techniques. The observed lines were taken 
from the following references, listed by spectral band with the aperture used (where given): UV 
lines from Kriss et al ( [1992b| ) (18"0); optical lines from Osterbrock & Parker ( |1965[ ) (trailing 
long slit), Anderson (|l970|) (8" x 8"), Wampler (jT97lD (10" 0), Koski ( |l978| ) (2.7" x 3.4") and 
Neugebauer et al ( |l980| ) (~ 10" x 20"); optical to near-IR lines from Shields & Oke ( g97p 
(10" x 10"), near-IR lines from Oliva & Moorwood ( |l990| ) (6" x 6"), Marconi et al dl996|) (long slit 



of width 4.4"), Osterbrock k Fulbright (|l99|) (3" x 18") and Osterbrock, Tran & Veilleux (|l99|) 



(long slit of width 1.2"); IR lines from Thompson ( |1996; ) (2" x 10") and the ISO-SWS (14" x 20" 
to 20" x 33"). 

As is discussed in detail by Alexander et al. ( |1999| ), the result of such a compilation is 
generally not self-consistent, and only a small subset of the of the ~ 120 lines can be used. First, 
we excluded low and medium excitation lines {E- lon <; 100 eV) observed through apertures whose 
smaller dimension is less than 3", so as to avoid significant loss of light due to incomplete coverage 
of the line emitting region. The observed emission from the high excitation lines is centrally 



concentrated in the inner < 4" (Marconi et al 1996 ), which are covered even by the smallest 



apertures used. Second, we excluded lines with E\ oa < 30 eV to avoid using lines that may be 
photoionized primarily by young hot stars or other non-AGN, lower-energy excitation processes. 
This also reduces the loss of light bias. Third, we excluded narrow lines whose measured flux 
is uncertain, either because the flux is very low (flux less than 5 x 10~ 13 erg s _1 cm -2 , ~2% of 
the strongest line), or because it has a significant broad component, such as the NvA1240 and 
CivA1549 lines. Fourth, we excluded the [Fex] A6734 and [Fexi] A7892 lines, whose collision 
strengths are highly uncertain and therefore cannot be modeled reliably. The final, much reduced 
line list includes 22 lines (Table |l|), which we use for obtaining the best-fit SED. Whenever more 
than one measurement of the line exists, we quote the average flux and use the rms scatter as 
an error estimate. Important IR lines that were not included in the final list were nevertheless 
compared to the best fit model predictions to verify that there are no gross inconsistencies. 



4. The photoionization models 

The method of constructing photoionization models for the NLR and, in particular, the 
"log 2 S fit procedure" for obtaining the best fitting SED and gas parameters is described in 
detail in Alexander et al. ( |1999| ), and is summarized here briefly. We parameterize the SED as a 
piece-wise broken power-law (Fig. [I]), and enumerate on the different possibilities of connecting 
the power-law segments. We test a large number of simplified NLR gas models, which consist of 



- 8 - 



optically thick (radiation bounded) clouds whose ionized surfaces partially cover a spherical shell 
around the continuum source. Each cloud extends in the radial direction as far as it takes to 
effectively absorb all the ionizing UV photons (specifically, until the hydrogen ionization fraction 
falls below 2%). The gas clouds are parameterized by their chemical composition, the hydrogen 
density n, the ionization parameter U at the irradiated face of the cloud, and the filling factor 
F. In addition, an asymmetry parameter A expresses the ratio between the ionizing flux directed 
towards the NLR and that directed towards the observer. This describes situations where the 
continuum source is not isotropic, or where only a fraction / re fl of the continuum is reflected 
towards the observer (A = l// re fl)- 

For each NLR gas model, the fit procedure uses the observed line fluxes to derive the best-fit 
SED for that gas model, and from it to derive in a self-consistent way the corresponding covering 
factor C, the inner NLR angular radius 6, the width of the ionized region A#i on (defined here as 
the radial extent of the Balmer lines emitting gas) and the reddening coefficient -E?b-v- These 
parameters are constrained by the observations and so can be used to limit the range of acceptable 
NLR gas models. The log 2 S fit procedure assigns a score S to the best-fit model, which means 
that the model line fluxes fit the observed ones up to a factor S, on average. The worst-fitting 
line and the factor by which it deviates from the observed value are also recorded. Monte-Carlo 
simulations are used to calculate confidence limits on the best-fit SED. In addition, we calculate the 
correlation between the model-to-data line ratios and the lines wavelengths, ionization potentials, 
critical densities and "deplitivity" (the tendency of an element to be depleted into dust). These 
residual correlations test whether the remaining inconsistencies in the best-fit model are related 
to inaccurate modeling of the reddening, the spectral hardness / softness of the continuum, the 
gas density or its dust content. The correct model should not display any such correlations. The 
final, global best-fit SED is the one with the best S'-score among all the NLR gas models that are 
consistent with the observed NLR geometry and reddening, whose worst-fitting line is not too far 
from the observed value, and which display no significant residual correlations. 

We assume in all models that the gas clouds have constant density and that / re fi = 0.005. We 
investigate two classes of models. The first consists of models with a single type of cloud. We 
enumerate on different values for the ionization parameter (log U = —1,-2 and —3), gas density 
in = 2000 and 10 4 cm -3 ) and filling factor (logF = —2,-3 and —4 at the ionized surface). We 
test two different sets of non-solar abundances, the low oxygen set and high nitrogen set (§ |2,2| ). 
We test two different possibilities for the radial run of the filling factor. The first is that F is 
constant, which corresponds to either a static distribution of clouds, a rotating distribution of 
clouds, a linear constant velocity outflow, a strongly decelerating outflow at a constant opening 
angle, or an outflow where clouds are continuously added to the flow (e.g. from the molecular 
torus). The second is F oc r -2 , which corresponds to a constant velocity and constant opening 
angle outflow where the clouds are formed at the base of the flow. The second class of models have 
two types of gas clouds, and are constructed by combining all the possible pairs of one component 
models. We assume that the clouds do not obscure one another. 
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The photoionization calculations were carried out using the numerical photoionization code 
ION99, the 1999 version of the code ION described by Netzer ( |1996j ). 

5. Results 

5.1. Single component models 

We find that single component models generally fail to fit the observed line ratios and the 
observed constraints on the geometry of the NLR. Low U models (log U = —3) cannot reproduce 
the high excitation lines regardless of the values of the other model parameters. For the F oc r -2 
models, the low U models significantly over-estimate the size of the NLR, while the high-C/ models 
require unphysical covering factors exceeding unity. Single component models with a high U and 
low constant filling factor do somewhat better, although such models are problematic because a 
constant filling factor is inconsistent with simple scenarios of NLR outflow. The best fit single 
component model has low oxygen abundance with \ogU = — 1, logn = 3.3 and log-F = —3. This 
model can reproduce the observed lines up to a factor of 2 and predicts 6 = 0.8", A0i on = 4.8", 
C = 0.35 and -Eb-v = 0.18. These values are roughly consistent with the observed constraints. 
However, this model under-predicts the observed OviA1035 flux by a factor of ~5 and shows a 
residual negative correlation with Ao at the 5% confidence level. Like all the models investigated 
here, the best fit SED displays a deep trough at 4 Ryd (log / = -27.4, -29.0, -27.4, -28.2 at 
2, 4, 8 and 16 Ryd, respectively). A fit of similar quality is obtained with the high nitrogen 
abundance set. 

5.2. Two component models 

A better fit to the observations is provided by the best-fit two component model (Table ||). 
This model fits the 22 observed line fluxes to within a factor of 1.9 on average (experience 
shows this is as well as one can expect for AGN photoionization models). The worst fitting 
line, [Arvi] A4.5/im, is under-predicted by a factor of 4. The model-to-data ratios of individual 
lines are displayed in Fig. ||. The low excitation IR lines {E lon < 30 eV) such as [Neil] A12.8^m, 
[S in] A18.7/UUL and [Sm] A33.5;um, which were not used in the fit, are nevertheless all consistent 
with the observations to within a factor of 2. However, the model-to-data line ratios for these lines 
are not much smaller than 1, as would be expected if there is a significant contamination from gas 
photoionized by starbursts. The high excitation IR lines such as [Six] A1.25yum, [Six] A1.4/um and 
[Siix] A2.6/im, which were not used in the fit, are also consistent with the observations to within 
a factor of 3. There are no statistically significant residual correlations between these ratios and 
Ao, E[ on , the deplitivity or n c . The best fit values of the extinction (2?b-V ~ 0.2) is in agreement 
with other estimates (§^]^). The covering factors are somewhat larger than is indicated by the 
observed opening angle of the emission cone (§ [2.ip . The dense compact component requires a 
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covering factor of C = 0.45, which corresponds to a bi-cone with an opening angle of ~115°. The 
less dense, more extended component requires a covering factor C = 0.26, which corresponds to 
a single cone with an opening angle of ~130° or a bi-cone with opening angle of ~90°. Figure ^ 
shows the best fit SED for this gas model together with the 99.9% confidence limits on it (the 
lower confidence limits at 2, 4 and 16 Ryd were not calculated as they extend beyond the SED 
grid). This SED shows the generic trough that is seen in the best-fit SED of all the models we 
investigated. This model has a less than solar oxygen abundance (Netzer |1997 ; Netzer &i Turner 



1997 ). We find only a few two-component models with high nitrogen abundance that come close 



to a reasonable fit to the observations. Of these, almost all display a negative correlation with 
deplitivity at the 5% to 10% significance, which may indicate that the line emitting gas is more 
depleted than is assumed by the high nitrogen abundance set. 

We conclude that best-fit two component model provides a better, but not an overwhelmingly 
better fit to the observed line ratios than the single component model. Its marked advantage 
over the single component models lies in its consistency with the observed NLR geometry and 
kinematics. This is further discussed in the next section. 



6. Discussion 



As was discussed in detail by Alexander et al. (1999), there are various degeneracies between 
the parameters that describe the gas model (n, A, F, U, C, 6, and A#i on ). These degeneracies 
allow the fit procedure to converge to a robust best-fit SED even when the assumed values of n, 
A, F or U significantly differ from their true values, since this can be compensated to a large 
degree by a suitable modification of the gas geometry (C, 0, and A0i on ). This property of the 
fit procedure is especially important in the analysis of NGC 1068, where observations indicate 
that the actual properties of the NLR gas are much more complex than can be modeled by our 
family of simplified gas models. For this reason, we place more weight on the fact that all the 
best-fit SED models, whether one or two-component, display a deep trough at 4 Ryd than on the 
determination of the exact values of the gas parameters. 

The trough in the SED is required for reproducing the relative line fluxes of the high and low 
ionization species. To check the robustness of this result, we attempted to re-fit the observed line 
fluxes with all of our one and two-component models using an approximate single power-law SED 
(log / = -25.8, -26.6, -27.4, -28.2 at 2, 4, 8 and 16 Ryd, respectively), which was held fixed in 
the fit procedure. In all cases the low excitation lines were over-estimated with respect to the high 
excitation lines, regardless of the values of U, n or F. For example, when the power-law SED is 
applied to the best-fit two-component gas model (Table |2|), the low excitation lines (Ei on <^ 50 eV) 
are over-estimated by the model by up to a factor of 13, whereas the high excitation lines 
(-^ion k; 100 eV) are under-estimated by up to a factor of 23. The overall mismatch of this SED 
with the observed line fluxes is reflected in both the poor fit score of S = 4.1 and in the very 
strong residual anti-correlation between the line ratios and -Eion; whose random probability is 
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10~ 4 . We caution against generalizing this result to mean that every AGN that exhibits a hard 
emission line spectrum has an absorbed ionizing SED. The emission line spectrum reflects the gas 
parameters, such as U, n and F, no less than it does the ionizing SED. It is necessary to have some 
knowledge of the likely range of values for these parameters in order to interpret the hardness of 
the line spectrum. Alexander et al. ( |1999| ) provide a counter-example where subtle cancellations 
between the SED and the gas properties lead to a situation where an AGN (NGC4151) with a 
hard absorbed SED has a softer emission line spectrum than another AGN (Circinus) with a soft 
unabsorbed Big Blue Bump. 

Although we do not claim to fix the gas parameters with certainty, the best fit model (Table |2|) 
is reassuringly consistent with the observations, which broadly indicate that the integrated NLR 
emission originates in two components. Component A of the model can be interpreted as a system 
of dense (n = 10 cm -3 ), centrally concentrated (0.3" < 9 < 1.4") outflowing gas clouds (F oc r~ 2 ) 
with a relatively high filling factor (F ~ 0.01) and high ionization parameter (logU = —1). 
Component B of the model can be interpreted as a more extended distribution (1.9" < 6 < 4.5") 
of lower density gas (n = 2 x 10 3 cm -3 ) with no net outflow (F = const.), with a lower filling 
factor (F = 0.001) and a lower ionization parameter (log U = —2). Component A contributes 
58% of the total line flux in 22 lines listed in Table |], with the contribution to individual lines 
ranging from 45% of the low excitation [Siv] A10.5^m (-Ei on = 34.8 eV) to more than 99.9% of the 
very high excitation line [Siix] A3.9/mi (Ei on = 303.2 eV). Its large covering factor indicates that 
there is probably a significant contribution of flux from the inner ~ 1" of the diffuse SW emission 
cone as well as from the bright NE cone. Component B contributes the remaining 42% of the 
total line flux, mainly in the lower excitation lines. Its covering factor is small enough for it to be 
concentrated mostly in the NE bright emission cone, as is observed. 

The best-fit procedure indicates that models with the low oxygen abundance set fit the 
observed line fluxes somewhat better than models with the high nitrogen abundance set. In 
particular, the Oiii] A1663, whose unusual relative weakness was an important argument for 
assuming non-solar abundances (Netzer |1997 ; Kraemer et al. 1998| ), is well reproduced by the 



best-fit low oxygen model with a model-to-data line flux ratio of 1.4 (Fig. |^), even though it was 
not included in our fit since it didn't pass the minimal flux criterion. 

The trough in the best fit SED (Fig. |2|) can be interpreted as an absorption trough due to an 
absorber between the continuum source and the NLR. The energy resolution of the SED template, 
which is limited by the computational cost of enumerating over all the SED combinations, is too 
low to allow detailed modeling of the absorber. Figure || shows an example of how a quasi-thermal 
big blue bump that is absorbed by neutral hydrogen would appear in our low resolution SED 
reconstruction. We find that the trough is consistent, for example, with an absorber that shadows 
the entire NLR (C a b s = 1) and has a column density of N H o = 6 x 10 19 cm -2 in neutral hydrogen, 
or with an absorber that allows a small leakage of unfiltered radiation (C a b s = 0.999) and a column 
density of N H o = 10 20 cm . A similar trough, consistent with an absorber of C a b s > 0.99 and 
N H o = 5 x 10 19 cnr 2 , was discovered in the reconstructed SED of Seyfert 1 galaxy NGC 4151 
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( Alexander et al. 1999). That absorber was also detected in the HUT absorption line spectra 



of the UV continuum of NGC 4151 (Kriss et al. [Tgggg , [Lggp . The HUT spectra of NGC 1068 



(Kriss et al 1992b| ) do not have a high enough S /N to allow the detection of absorption lines in the 



scattered UV continuum of this AGN or against the stellar background (G. Kriss, private comm.). 
We predict that future sensitive absorption line studies should reveal the presence of such an 
absorber. 

The bias in our results due to the fact that we neglected the line emission from the absorbing 
gas is likely to be small if the absorber is similar to the dense, high velocity UV absorber that was 
detected in NGC 4151. Such an absorber will not emit forbidden lines, and its permitted lines 
will be broader than typical NLR lines. Only 3 of the 22 lines we used in our fit are permitted 
lines, and we did not use lines that are contaminated by broad components. A highly ionized and 
optically thin absorber will produce strong OviA1035 line emission in excess of the typical NLR 
emission. It is therefore interesting that unlike the two forbidden [Om] A5007 and [Oiv] A25.9/um 
lines and the semi-forbidden Om] A1663 line, which are well reproduced by the best fit model, the 
observed OviA1035 line is 3.2 stronger than predicted (Fig. |3|). 

We have, up to now, applied our SED reconstruction method to ISO-SWS observations of IR 



coronal lines of three AGN: the Seyfert 2 Circinus galaxy (Moorwood et al. 1996; Alexander et al 



I1999D , the Seyfert 1 galaxy NGC 4151 (Alexander et al. |1999D , and the Seyfert 2 galaxy NGC 1068 



(this work). In one of these (Circinus), we detect a Big Blue Bump that peaks at p> 50 eV. In the 
other two we detect deep troughs, which are consistent (but not exclusively so) with a Big Blue 
Bump that is absorbed by neutral gas interior to the NLR. Our findings thus far are consistent 
with the picture that luminous Seyfert galaxies are powered by thin accretion disks that produce 
a quasi-thermal Big Blue Bump, and that in a large fraction of them the NLR sees a partially 
absorbed ionizing continuum, as suggested by Kraemer et al. ( |1998 ). 
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Fig. 1. — The template used for enumerating on the SED of NGC 1068. The most and least luminous 
SEDs are indicated by the lines, together with the relative sense of their spectral hardness between 
1 and 150 Ryd. 

Fig. 2. — The best fit SED of the two-component model (Table ||). 

Fig. 3. — The model to data line ratios for the best fit two-component model (Table §). The 
O in] A1663 line ratio (white circle), which was not used in the fitting procedure is also displayed. 

Fig. 4. — An example of how a quasi-thermal Big Blue Bump that is absorbed by neutral hydrogen 
would appear in the low resolution reconstructed SED. The unabsorbed bump (thin line) and the 
bump after absorption by a C a b s = 1, N H o = 6 x 10 19 cnr 2 absorber (long dashed line) and a 
Cabs = 0.999, N H o = 10 20 cm -2 absorber (short dashed line) are superimposed on the best fit SED 
model. 
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Table 1. The compiled emission line flux list. 



Line A £ ion a f t h A/ £ c 

/im eV 10 -13 ergs _1 cm -2 



UV, optical and NIR lines 



Ovi 


0.1032+1038 


113.9 


37.4 


3.1 


Niv] 


0.1487 


47.4 


5.1 


1.1 


He II 


0.1640 


54.4 


17.7 


5.2 


[Ne v] 


0.3426 


97.1 


15.7 


6.7 


[Ne in] 


0.3868+3969 


41.0 


19.2 


2.6 


He II 


0.4686 


54.4 


6.15 


1.50 


[Om] 


0.4959+5007 


35.1 


256. 


27. 


[Si vi] 


1.96 


166.8 


7.5 


0.5 


[Si vii] 


2.48 


205.1 


8.3 


0.9 


ISO-SWS IR lines 


[Mg vm] 


3.028 


224.9 


11. 


1.1 


[Si ix] 


3.936 


303.2 


5.0 


0.6 


[Mgiv] 


4.49 


80.1 


7.6 


1.5 


[Ar vi] 


4.528 


75.0 


15. 


3. 


[Mg vii] 


5.50 


186.5 


13. 


? 


[Mgv] 


5.61 


109.2 


18. 


2. 


[Ne vi] 


7.64 


126.2 


110. 


11. 


[Siv] 


10.51 


34.8 


58. 


6. 


[Nev] 


14.32 


97.1 


97.0 


9.7 


[Nem] 


15.56 


41.0 


160. 


32. 


[Nev] 


24.32 


97.1 


70. 


7. 


[Oiv] 


25.89 


54.9 


190. 


20. 


[Ne in] 


36.01 


41.0 


17. 


3. 



a The ionization energy required to produce the 
emitting ion from the preceding ionization stage. 

b The observed flux. For permitted lines, the flux of 
the decomposed narrow component is quoted. 

c The error estimate on the observed flux. The ISO- 
SWS errors are estimated by the scatter in various 
methods for defining the underlying continuum and 
measuring the line. The errors do not include systematic 
calibration errors, which are generally smaller than 
30%. A question mark means that error estimates are 
unavailable. 



Table 2. The best-fit two-component model. 



Component A Component B 



Input parameters: 



Best-fit results: b 



Residual correlations: 



U 
n 

Composition 
A 

pa. 



i log r 



s 

max Si c 
worst line 
log/ 2 d 
log/4 

log/8 

l0g/l6 

C 

e 

A#ion 

Eb-v 

Qion 

l0gQion/s _ 

(h V y 

A 

p, 

depletion 



1 o 



0.1 

cm" 
low oxygen 
200 
I0~ 2 h 

-2 



0.01 
10 3 ci 
low oxygen 
200 
10- 3 h 




1.9 

4.0" 1 
[Arvi] X4.5^m 
-29.0 (-29.0, -26.6) 
-29.0 (-29.0, -26.6) 
-27.4 (-27.4, -26.6) 
-29.0 (-29.0, -28.2) 
0.45 0.29 
0.26" 1.9" 
1.1" 2.6" 
0.23 (0.18, 0.29) 
0.24 s _1 cm- 2 
54.2 
5.1 Ryd 



Correlation 
-0.19 
0.00 
+0.11 
+0.03 



Random prob. 
0.23 
1.00 
0.70 
0.86 



a Assuming h = 0.65. 

b Values in parentheses are the 99.9% confidence intervals. 
c Model to data ratio, given as the reciprocal when < 1. 
d Flux in erg s 1 cm 2 Hz 1 . 

e Ionizing photon luminosity assuming isotopic emission. 
f Mcan ionizing photon energy between 1 and 16 Ryd. 
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